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Abstract. We conduct 3D magneto-hydrodynamic (MHD) simulations in order to test the 
Oh' stability of the magnetic equilibrium configuration described by Duez & Mathis (2010). This 
analytically-derived configuration describes the lowest energy state for a given helicity in a stellar 
radiation zone. The necessity of taking into account the non force-free property of the large- 
scale, global field is here emphasized. We then show that this configuration is stable. It therefore 
provides a useful model to initialize the magnetic topology in upcoming MHD simulations and 
stellar evolution codes taking into account magneto-rotational transport processes. 
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1. Introduction 

The large-scale, ordered nature of magnetic fields detected at the surface of some Ap, 
and B type stars and the scaling of their strengths (according to the flux conservation 
scenario) favour a fossil hypothesis, whose origin is not yet elucidated. To have survived 
since the stellar formation, a field must be stable on a dynamic (Alfven) timescale. It was 
suggested by Prendergast (1956) that a stellar magnetic field in stable axisymmetric equi- 
y—i • librium must contain both poloidal (meridional) and toroidal (azimuthal) components, 
since both are unstable on their own (Tayler 1973; Wright 1973). This was confirmed 
recently by numerical simulations (Braithwaite & Spruit 2004; Braithwaite & Nordlund 
1^ | 2006) showing that an arbitrary initial field evolves on an Alfven timescale into a sta- 
■ ble configuration; axisymmetric mixed poloidal-toroidal fields were found. On the other 
hand, magnetic equilibria models displaying similar properties have been re-examined 
analytically by Duez & Mathis (2010). We here address the question of their stability. 



2. The model 



We deal with non force-free magnetic configurations (i.e. with a non-zero Lorentz 
• »"H . force) in equilibrium inside a conductive fluid in absence of convection. Several reasons 
inclined us to focus on non force-free equilibria. First, Reisenegger (2009) reminds us 
that no configuration can be force-free everywhere. Although there do exist "force-free" 
configurations, they must be confined by some region or boundary layer with non-zero or 
singular Lorentz force. Second, non force-free equilibria have been identified in plasma 
physics as the result of relaxation (self-organization process involving magnetic reconnec- 
tions in resistive MHD), e.g. by Montgomery & Phillips (1988). Third, as shown by Duez 
& Mathis (2010), this family of equilibria is a generalization of Taylor states (force-free 
relaxed equilibria; see Taylor 1974) in a stellar context, where the stratification of the 
medium plays a crucial role. The equilibrium obtained is described in detail in Duez & 
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Mathis (2010) as the lowest energy configuration conserving the invariants of the prob- 
lem (during the relaxation phase) which are the magnetic helicity (preventing the rapid 
energy decay) and the mass enclosed in magnetic poloidal flux surfaces (to account for 
the non force- free property). 

3. Stability numerical analysis 

We use the Stagger code (Nordlund & Galsgaard 1995). We model the star as a self- 
gravitating ball of ideal gas (7 = 5/3) with radial density and pressure profiles initially 
obeying the polytropic relation P oc ; with index n = 3, therefore stably stratified. 

More details on the numerical model setup can be found in Braithwaite & Nordlund 
2006. The configuration is then submitted to white perturbations (1% in density). The 
dynamical evolution of the mixed configuration is compared to its purely poloidal and its 
purely toroidal components, whose behaviour are well known to be unstable due to kink- 
type instabilities. The magnetic and velocity amplitudes are plotted on Fig. 1. As we can 




Figure 1. Time evolution of the (log) amplitudes in azimuthal modes m = to 4 averaged over 
the stellar volume of the magnetic field (top row) and the velocity field (bottom row) in the 
simulations with the purely poloidal field (left), purely toroidal field (middle) and the mixed 
field (right). Initially, all the magnetic energy is in the m = mode since the initial conditions 
are axisymmetric. 

see, in contrast to these unstable configurations, the mixed poloidal-toroidal one does not 
exhibit any sign of instability, even for high azimuthal wavenumbers (up to about 40). 
This is the first time the stability of an analytically-derived stellar magnetic equilibrium 
has been confirmed numerically. This result has strong astrophysical implications: the 
confuiguration provides a good initial condition to magnetohydrodynamic simulations 
and magneto-rotational transport to be included in next generation stellar evolution 
codes; furthermore it could help to appreciate the internal magnetic structure of neutron 
stars and magnetic white dwarfs. 
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